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In this paper, we report our first attempt at determining NN potentials in the parity odd sector 
including the spin-orbit force in lattice QCD, employing the method to extract successfully parity 
even NN potentials from Nambu-Bethe-Salpeter (NBS) wave functions through the Schrodinger 
equation. Using Nf — 2 CP-PACS gauge configurations on a 16 3 x 32 lattice at a ~ 0.16 fm 
and nij ~ 1 . 1 GeV, we calculate central, tensor and spin-orbit potentials in the parity odd sector. 
Although statistical errors are still large, we observe that the qualitative features of these potentials 
roughly agree with those of phenomenological potentials. 



X 



The XXIX International Symposium on Lattice Field Theory - Lattice 2011 
July 10-16, 2011 

Squaw Valley, Lake Tahoe, California 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Nuclear forces in the odd parity sector and the LS forces 



Keiko Murano 



1. Introduction 

The nucleon-nucleon (NN) potential is one of the fundamental quantity to study various prop- 
erties of atomic nuclei. In the past few decades, based on various models[|I], ^, ^ or the chiral 
effective field theory several realistic NN potentials have been constructed to reproduce the 
NN scattering phases of more than 4000 data points with % 2 //Vdof ~ 1- On the other hand, quite 
recently, a new method to extract NN potential in QCD has been proposed[Q]. For the moment, 
the method has been successfully applied to the potentials in the parity even sector at the leading 
order of the derivative expansion, i.e., the central force ( V , V<j ) and the tensor force (Vy ) 
are obtained for various cases including NN, YN, and YY systems^, 0, [sj] (where Y represents a 
hyperon) and meson-baryon system^]]. In addition, the method is recently extended to investigate 
3N forces [|l0|]. In contrast, the method has not yet applied to the forces in the parity odd sector as 
well as the spin-orbital force. They are needed for the complete determination of NN potentials 
at next-to-leading order (NLO). In particular, the spin-orbit force, even though it appears at NLO 



of a derivative expansion, is important, as is indicated by an analysis of phase shift|ll|, [12p. It 
is also important in explaining the ls-splitting of the (hyper) nuclear spectra and the magic num- 
bers in nuclei, i.e., it induces the one-body spin-orbit term in the average single-particle nuclear 
potentially, 



The purpose of this paper is to extract potentials in the parity odd sector including the spin- 
orbit force for the NN system. For this purpose, we construct NBS wave functions which include 
higher angular momenta, in order to determine the central force, the tensor force and the spin-orbit 
force in the parity odd sector (Vq , V„ \ Vj \ V^ s ^). Our calculation has been made by using 
Nf = 2 CP-PACS configurations on a T x L 3 = 16 3 x 32 lattice at j8 = 1.95 (a ~ 0.16 fm) and 
K = 0.1375 (m n ~ 1100 MeV). 

2. Extraction of potentials 

2.1 Definition of potentials 

In the method, the NN potential is constructed form Nambu-Bethe-Salpeter (NBS) wave func- 
tion in the center of mass (CM) frame, defined by 

<P a ^r;E) = (0\p a (x)n p (y)\p(+k)n(-k)}, (r = x-y) (2.1) 

where p a and na denote local composite nucleon operators with spinor indices a, /3, and E is 



kinetic energy related to the relative momentum k = \k\ as E = 2 y k 2 + mjf — 2m;v(~ k 2 /m^) with 
rriN being the nucleon mass. We define NN potentials from the NBS wave function below the 



inelastic threshold through the following Schrodinger type equation \ \15[ |16| ] 

V 2 \ r i 

E + — )(j)(r,E)= P + v( +) (r)+p-v(->(r) <j>(r;E), (2.2) 
m N J L J 

where P + (P~) denotes a projection operator for parity even (odd) with the corresponding potential 
y( + ), (V^), which is decomposed as 



V^(r) 



V ( r ) + V e ( r ^i • & + V T ( r )Si2 + V£> (r)L -S+ (NNLO) 



(2.3) 



2 



Nuclear forces in the odd parity sector and the LS forces 



Keiko Murano 



(±) 



r(±) 



with the central and spin-dependent central forces V and , the tensor force Vj , and the spin- 
orbit force V^y . While V , vjj and Vy are of leading order(LO) in the derivative expansion 
of non-local potentials, appears at next-to-leading order(NLO). Once the above NBS wave 
functions are calculated in lattice QCD simulations, potentials can be extracted by solving eq.(2.2). 

2.2 Spin-singlet potentials 



For the spin-singlet and parity odd sector, the Schrodinger type equation eq.(2.2) is reduces to 



v 2 \ .(-) 



E + — ) Co(^) = Vc { s J (r)<l>^(r,E), 
m N l 



(2.4) 



where tp^J (r;E) = p( ^P s ~°(j)(r,E) with the spin-projection operator P s and Vc s J (r) = V (r) — 
(r) denotes the central potential in the spin-singlet channel. From this equation, we then obtain 



Vc [ s J (r)=E + 



1 V 2 0(r;£) 



m N (j>(r;E) 

2.3 Spin-triplet potentials including spin-orbit force 



(2.5) 



For the spin-triplet and parity odd sector, the Schrodinger type equation eq.( |2.2| ) reads 



v 2 \ 



m N J L J 



(2.6) 



Where ^ = {(r) = />H P 5 " 1 <p (r; E) and Vc^J^r) = V (r) + V~(r). Note that the isospin of these 
V ~ and Vg is different from that of the spin singlet sector. 

In order to determine the spin-orbit potential for parity odd sector, we consider three linear 
independent set of NBS wave functions, <p Xi with i = 1,2,3, constructed by 

( (, x -(r)=P Xi P s=1 P^^(r;E). (2.7) 

where P x < denotes a projection operator to state X,-, labeled by the orbital angular momentum L and 
the total spin /. Eq. (2i) is then decomposed into three independent equations as 

V 2 \ 



E + —\ <b Xl (r)=P x > 
m N J 



VP (r) + Vf>(r)S n + V^>(r) L S <j> (r) 



,(-) 



(2.8) 



By solving these equations, potentials Vq ,Vj and vj^ ^ can be obtained as 



V 2 /m N 0* 2 (r) 

72 / m „ AX}(7!\ 



/Vt ) (r)-E\ 


1 / 


vjr\r) 


= M(r)- 1 




V 



(2.9) 



where Mir) denotes 3x3 matrix defined by 

M(f) 



/ 0* (?) P x ^S l2 $ (r) P x> L ■ S(j) (r) \ 

^-(r) P X2 S n ${r) P Xl L-S§{r) 
^ Xi {r) P Xi S n ${r) P Xi L-S§(r) j 



(2.10) 
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3. Construction of NBS wave functions in lattice QCD 

The NBS wave function is obtained from a 4-point nucleon correlation function on the lattice 



as 



1 



G a ,p(x-y,t-t ; a f) = jj J £(°\ T Pa(x + 7,t)n p (y + r,t)^ J - J ' S - s (t ) |0) (3.1) 



K;p J,S ~ S (x-y> E °) (p(+ko)»(-ko) I J 



J=J~S=S 



(0)\0)e 



-E (t-to) 



t-tQ > 1, 



(3.2) 



where the summation over r is performed to select the two nucleon system with vanishing total 
spatial momentum. For p(x) and n(y), we employ the following local composite operators 



p(x) = e a bc{u' a {x)Y 5 db(x))u c (x), n(x) = e afoc (M^(x)Cy 5 ^(x))J c (x), 



(3.3) 



where a, b and c denote color indices, while ^f J=J ' s=s denotes a two-nucleon source operator with 
a definite total spin S = S(S = 0, 1) and a definite total "angular momentum" / = /, which is indeed 
an irreducible representation of the cubic group as f = A\ ,A2,E, T\ , T 2 . 

For such two-nucleon sources, we take two-nucleon momentum wall sources defined by 



where 



Paif) 
Np{f) 



Jap(f)=Pa(f)Np(f), 
£ £abc (" fl (xi)Cy 5 <4(x 2 )) M C)Ce (x 3 )/(x 3 ] 

xi,x 2 ,x 3 

£ £abc {u a {xi)Cy 5 d h {x 2 )) d Cjj} (x 3 )f(x 3 ) 



(3.4) 



(3.5) 



Xl,X 2 .X 3 



with / being one of the following source functions, each of which corresponds to a plain wave 
parallel to one of spatial coordinate axes as 



/(°)(r) = e X p(-27nx/L), f l >(r) = exp(-2^y/L), f z >(r) = exp(-2^z/L), 
/ 3 '(r) =exp(+2nix/L), f w (r) = exp(+2niy/L), / 5 '(r) = exp(+2niz/L). (3.6) 

Note that an element g in the cubic group O with 24 elements acts on these plane waves as 

f^^Uijig)^, (3.7) 
j 

where the 6 x 6 representation matrix U (g) are generated by representation matrices for c y and c z , 
which are the rotations by 90 degree around y and z axes, respectively, and are explicitly given by 



U(Cy) = 



"0 1 0" 




"0 1 0" 


10 




10 


10 
1 


, U(c z ) = 


10 
10 


10 




10 


.1 0. 




.0 1. 



(3.8) 
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Note also that the spatial reflection for these source functions is represented by the complex conju- 
gation. By a cubic group analysis, the orbital part of this momentum wall source is decomposed into 
A\ © E + © . Therefore, for the parity odd sector, J p we can access is (L = T{~ ) ® (S = A i ) = Tf 
for the spin-singlet sector and (L = Tf ) ® (S = T\ ) = A j~ © E ~~ © Tf © T 2 ~ for the spin-triplet sector. 
The wall source with the definite total angular momentum is now constructed as 

Afv®) =ll x^nM/^if^^aiB- 1 )^^ 1 ) (3-9) 

z ^ g£0 

where and X is) denote the dimension and the character for the irreducible representation 
J, respectively. Hereafter the Dirac indices a,j8 are restricted to upper components (in the Dirac 
representation). The total spin S is projected by the spin projection operator as 

p(s=o) = (l_ ai -c72)/4 and p( s=1 ) = (3 + 0\ ■ c-ijjA. Finally, the parity projection is defined by 

P (±) ■ Safitf®) = \ ( JapV®) ± A/3(/ 0> )) • (3-1D 

Although the orbital part for our two-nucleon source operators in the parity odd system consists 
of L = T\ only, other orbital components of the NBS wave function at the sink are induced through 
the effect of the tensor force in the case of the spin-triplet sector. In this paper, for simplicity, we 
consider only an L = T\ component of the NBS wave function as 

(j> J L Z J T f= S (r;E ) = £ X iTl) (g)*<l> J=f > S = § (R(g)r;Eo). (3.12) 
geo 

where R(g) denotes the rotation matrix for a cubic group element g, which rotates only the orbital 
components of the NBS wave functions as x t- > x' = R(g)x. 



4. Lattice setup 

Our calculation is performed on a set of Nf = 2 dynamical QCD configurations generated by 
the CP-PACS Collaboration on a 16 3 x 32 latticejl7|], who employs the 0(a) -improved Wilson 
quark action with C$w = 1-53 at k uc i = 0.1375 and the RG improved gauge action (Iwasaki action) 
at j8 = 1.95. This set of parameters corresponds to the lattice spacing a = 0.1555 fm and pion 
mass m % = 1136 MeV( nucleon mass = 2165 MeV). The 4-point nucleon correlation function 
eq. ( |3~l| ) is calculated with the periodic and the Dirichlet boundary conditions along the spatial 
and the temporal directions, respectively. The spatial momentum of each nucleon is discretized as 
ki ~ 27irii/L(rij € Z) in the periodic BC. Since k = state is forbidden in the parity odd system, the 
ground state energy with parity odd system is not E = but E ~ (2n/L) /m^. In our analysis, we 



thus employ a free value E = (2k /L) /m^ = 1 15 MeV for an energy E in eq. (2.9). 
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Figure 1: (Left) The central potential in the spin-singlet and parity odd sector, calculated from P\ NBS wave 
function. Data at t — to = 7,8,9, 10 are simultaneously plotted. (Right)The central potentials in the spin- 
triplet and parity odd sector, calculated from 3 Po, 3 'Pi and ^Pi NBS wave functions. Data at t — to = 4,6,8 
are simultaneously plotted. 



5. Numerical results 

For the spin singlet sector, we calculate the J p = Tf NBS wave function with L = T\, ^[ZtI S_ °> 
whose dominant component corresponds to the l P\ NBS wave function. Using eq.(2.5), we obtain 



the central potential in this channel. Preliminary results at t — to = 7 — 10 are plotted in Fig. |T] 
(Left). We observe that the central potential has a strong repulsion at short distance (r < 0.5 fm) 
without attraction at longer distance. 

J— J S— 1 ~ 

For the spin triplet sectors, we calculate three NBS wave functions ((> L _f with / = A\,T\ 
and E, whose dominant components correspond to 3 Po, 3 Pi and 3 P2, respectively. Using these 
NBS wave functions, we extract spin-triplet potentials in parity odd sector. The central potential 
Vci_j (r) is shown in Fig. [T] (Right). Moreover the tensor potential Vf (r) and the spin-orbit force 
V L £ (r) are given in the left and right panels of Fig. ||, respectively. Large fluctuations at r ~ 1.2 
fm are caused by spatial boundaries at L/2 = 8, which breaks the rotational symmetry, while large 
statistical errors of the tensor force in the parity odd sector may be mainly caused by the fact that its 
magnitude is significantly smaller (few MeV) than other potentials (few-hundred MeV). Although 
systematic and statistical uncertainties are still large in our preliminary results, we observe the fol- 
lowing qualitative features. (1) Vc is repulsive at all distance with a repulsive core at short distance. 
(2) Vj is positive and quite small. (3) Vls is large and negative. These features qualitatively agree 
with those of phenomenological potentials^]. 



6. Conclusions 



In this paper, we have made a first attempt at determining NN potentials in the parity odd 
sector including the spin-orbit force. Using Nf = 2 CP-PACS gauge configurations on a 16 3 x 32 
lattice at a ~ 0.16 fm and m n ~ 1100 MeV, we have reported our preliminary results on central 
forces, the tensor force and the spin-orbit force in the parity odd sector. Although statistical errors 
are still rather large, we have observed that the qualitative behaviors of these potentials roughly 
agree with those of phenomenological potentials. 
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Figure 2: Potentials with spin-triplet state in parity odd sector, calculated from 3 Pq, i P[ and 3 P2 NBS wave 
functions. Data at t — to = 4,6,8 are simultaneously plotted. (Left) The tensor potential. (Right) The spin- 
orbit potential. 



Numerical calculations are performed on University of Tsukuba Supercomputer system (T2K). 
This work is supported by the Grant-in-Aid for Scientific Research on Innovative Areas(No.2004: 
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